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We present high resolution ARPES data of the bilayer superconductor Bi2Sr2CaC'u20s+s 
(Bi2212) showing a clear doubling of the near Ef bands. This splitting approaches zero along 
the (0,0) —> [iTjTt) nodal line and is not observed in single layer Bi2Sr2CuOe+5 (Bi2201), suggest- 
ing that the splitting is due to the long sought after bilayer splitting effect. The splitting has a 
magnitude of approximately 75meV near the middle of the zone, extrapolating to about lOOmeV 
near the (tt, 0) point. 



One of the central features of the cuprate superconduc- 
tors is that their physical properties depend strongly on 
the number of Cu — O planes (n) per unit cell. For exam- 
ple, increasing n from 1 to 3 almost universally results in 
large enhancements in the superconducting Tc [Q. This 
would seem to argue for the importance of some sort of 
coupling between the Cu — O planes within a unit cell. 
Despite this, there have been essentially no direct obser- 
vations of the intraccU coupling between the planes, and 
the vast majority of theoretical models for the cuprates 
focus on a single Cu — O plane, i.e. on the limit of zero 
intra and intercell c-axis coupling. 

In this letter, we present the first high en- 
ergy/momentum resolution angle resolved photoemis- 
sion (ARPES) data of a cuprate superconductor which 
shows a clear doubling of the bands in the normal 
state. The splitting is observed in the double layer com- 
pound Bi2Sr2CaCu20s+s (Bi2212) but is much weaker 
or absent in the single layer compound Bi2Sr2CuOe+s 
(Bi2201). The splitting approaches zero along the 
(0,0) (7r,7r) nodal fine and is maximal at M(tt,0), 
and appears to be larger for overdoped samples. These 
trends point towards an intracell c-axis coupling between 
the Cu — O planes (bilayer splitting) as the hkely origin of 
the splitting, although other effects such as microscopic 
phase separation may play a role. In addition to the new 
physics of this coupling, our new findings may help re- 
solve the recent controversy about the electron-like (^-^] 
and hole-like [Ij-Q Fermi surfaces which have recently 
been measured on Bi2212. 

The experiments were carried out at the Stanford Syn- 
chrotron Radiation Laboratory (SSRL), Stanford, CA, 
at the Advanced Light Source (ALS), Berkeley, and at 
the Synchrotron Radiation Center (SRC), Stoughton, 
WI. At each facility we used a Scienta 200mm energy 
analyzer running angle mode to simultaneously collect 
up to 84 individual Energy Distribution Curves (EDCs) 
along an fa 14° angular slice. The angle resolution was 



about 0.08° along the angular slice (the 9 direction) and 
about 0.4° in the perpendicular direction (the (j) direc- 
tion). At the experimental photon energy hv = 2A.leV, 
the Brillouin Zone edge (tt/o) will be k, 20° and the cor- 
responding momentum resolution will be (Afca;,Afcj^) w 
(0.0087r/a, 0.047r/a) with a the lattice constant. The en- 
ergy resolution was better than 20meV determined from 
the 10 — 90% width of a gold reference which was in elec- 
tric contact with the sample and was used to determine 
the Fermi Energy Ep. Data is shown from a lightly over- 
doped Bi2212 sample with a Tc = 85K and a lightly 
overdoped Bi2201 samples with Tc ~ 5K. The photon 
polarization direction was along the 9 direction and the 
r(0, 0) — Y{Tr, tt) high symmetry line. All data shown 
here was normalized by using the spectral weight in the 
high binding energy window {—O.SeV, —O.GeV) only H. 
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FIG. 1. (a),(b) ARPES intensity map within ±5meV of 
the Ef from the normal state (T — lOOK) of an overdoped 
Bi2212 sample. The photon energy was 24.7eV^ and the po- 
larization 40% along the 6 direction and 60% out of plane. An 
overlay of the data with the Bi2212 Brillouin Zone is shown in 
(a). The hole-like main FS (labeled "Main") plus superstruc- 
ture bands (labeled "S.S.") was obtained at 19eV by Ding et 
al. [5] (yellow lines). The FS crossings determined from MDC 
fitting are indicated by white dots. 



Figure la shows the schematic plot of the commonly 
accepted Fermi Surface (FS) topology of Bi2212 in 
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part of the first Brillouin Zone. The thick yellow lines 
are the main FS and the thin yellow lines are the FS 
produced by the crystal distortion (the superstructure, 
labeled "S.S."). Figure lb is the intensity distribution 
of the ARPES spectral weight at Ep in the normal state 
of a lightly overdoped Bi2212 sample, with the k space 
location of this panel indicated in figure la. Modulo ma- 
trix element effects, the high intensity locus at Ep should 
represent the FS topology In panel b, one can see 
that the main FS splits into two branches, one produced 
by band A and the other by B. This splitting behavior 
can also be seen in the superstructure FS, even though 
the intensity is much weaker there. This behavior is not 
expected within the generic FS topology shown by the 
yellow lines in figure la. 

In order to examine the splitting of these bands more 
carefully, we examine the /c-dependence of the states as a 
function of binding energy. Figure 2a shows data along 
the blue line {</) = 7°) in figure la, with the main band 
crossing Ep near the center of the window. The so-called 
"main band", which shows up as the high intensity fea- 
ture, disperses from binding energy « —0.12eV (around 
emission angle = —6°) towards Ep. Around 9 = —12°, 
another weaker feature appears at slightly deeper binding 
energy. This two-peak structure (labeled "A" and "B" 
here) can be easily seen in the energy distribution curve 
(EDC) shown in panel b, which was taken a± 6 = —12° 
(vertical white line in panel a) . We emphasize that such 
two-peak structure in EDC is self-evident without any 
analysis on the data. However, further quantitative de- 
termination of the peak positions, FS crossings, or split- 
ting is difficult from the EDCs, as a detailed understand- 
ing of their line shape and background is still unclear. In 
addition, the Fermi function distorts the line shape near 
the Fermi energy. 

On the other hand, it is now well established that for 
relatively dispersive states the Momentum Distribution 
Curves (MDCs, intensity versus emission angle 6 with 
fixed binding energy) have a simple line shape which can 
be modeled by a simple Gaussian or Lorentzian function 
pO|-p^. This is equivalent to saying that the electron 
self-energy is essentially fc-independent over the narrow 
/c-range of a MDC peak. A further advantage of the MDC 
method is that it won't be affected by the Fermi func- 
tion cutoff, so in principle one can track ARPES features 
all the way up to or even above Ep. Figure 2c shows 
the MDCs from the data of 2a over the full 9 range and 
from IQmeV above Ep to GOmeF below. The MDC at 
Ep is further shown in panel d (blue dots), along with 
a fitting (red line) of the main band into two separate 
Lorentzian peaks A and B and the superstructure band 
into two corresponding features (green lines), plus a lin- 
ear background (black dash line). For the fit, the angular 
width, position, and intensity of each feature was allowed 
to vary. It was found that the angular width of both 



branches were essentially identical, and the intensity ra- 
tio and the angular splitting between the branches found 
to be roughly constant as a function of binding energy 
(see figure 2e). The high quality fits allow an accurate 
determination of the peak positions for each component. 
Such fits were performed for each MDC at Ep between 
the angles of 3° and 8° [Q. In this way the white 
dots of figure la were determined which, as expected, 
perfectly match the high intensity locus. 




binding energy emission angle 6 x10' 

Figure 2 

FIG. 2. (a) False color plot oi E ys for the = 7° cut 
(blue line in figure la. (b) EDC at ^ = —12° from panel 
a (vertical white dash line). Two distinct features, A and 
B, can be clearly seen in this EDC. (c) MDCs from panel a 
from binding energies -|-10meV to — 60meV, with 5meV step, 
(d) The MDC at Ep (blue dots), including a deconvolution 
of the main band (red line shows the fitting result) into two 
Lorentzian functions A and B (green lines) plus two corre- 
sponding features in superstructure band (green lines) on top 
of a linear background (black dash line), (e) The energy de- 
pendence of the 6 value of MDC peaks A (close circles) and B 
(open circles). The error bar from the fitting is smaller than 
the symbol size. 

The energy splitting of these two bands can be most 
accurately found by fitting the MDC peaks for a range 
of energies. Figure 2e shows the results of the fits from 
the (j) = 7° data of 2c. The black dots are the peak cen- 
troids {9, in degree) for each of the components of the 
main band. The energy splitting can then be measured 
as the separation between the two MDC branches, which 
is immediately seen to be ^ 75mel^. 

There are several possible origins for the band split- 
ting, both intrinsic and extrinsic. We begin by discussing 
and ruling out the extrinsic effects. First among these is 
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the possibihty that the sample surface was slightly mul- 
tifaceted, which each facet contributing to a feature at a 
slightly different angle. This is not the case, as a post- 
cleave laser reflection plus optical inspection revealed an- 
gular deviations no greater than 0.2°, while the angular 
splitting of the states is on the order of 2°. X-ray Laue 
and LEED reflections showed sharp spots as well, indi- 
cating high quality crystalline of the samples. We also 
considered the possibility that there may be Bi2201 in- 
tergrowths in the Bi2212 samples, which would be very 
difhcult to detect from magnetization or x-ray measure- 
ments. In this possibility, it could be imagined that one of 
the branches was due to Bi2212 and the other to Bi2201. 
One could possibly imagine this effect causing the split- 
ting on one cleave, but not on the numerous cleaves we 
have studied. Also, we have cooled a sample down to 
35K which is below Tc for Bi2212 and above the Tc for 
Bi2201. Since the superconducting gap was observed in 
both branches, it rules out the possibility of the Bi2201 
contamination. 
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Figure 3 

FIG. 3. (a) MDC at Ef for the cl> = 7° cut of lightly over- 
doped Bi2201, under the identical experimental conditions as 
figure 2. (b) same data as figure lb except with the polar- 
ization 66% along the 9 direction and 34% out of plane, (c) 
MDC at Ef for the (f) = 7° cut (white line in figure 3b). 

Two main intrinsic possibilities exist: (a) the bilayer 
splitting due to the intracell c-axis coupling, and (b) mi- 
croscopic phase separation into hole rich and hole poor 
regions each with distinct bands and Fermi surfaces. We 
first discuss option (a), which we believe to be much more 
likely. We first remind the reader that in the bilayer 
splitting scenario, the electronic states from each of the 
two CuO planes per unit cell will couple, breaking their 
degeneracy. One set (antibonding) will increase in en- 
ergy while the other (bonding) will decrease. Within this 
picture we make the assignment on the figures as peak 
A=antibonding and peak B=bonding. Important infor- 
mation is garnered by comparing the spectra to that of 
the single layer compound Bi2201, which can not exhibit 
bilayer splitting. Figure 3a shows an MDC at Ep from 
the (j) = 7° angular slice of a lightly overdoped Bi2201 



sample taken under the identical experimental configu- 
ration as that of figure 2. Unlike figure 2d, a single 
Lorentzian function alone is enough to produce a very 
good fit to either the main or superstructure band. This 
indicates that there is no splitting of the bands in Bi2201, 
at least not one with a similar magnitude as observed 
in Bi2212. This Bi2201 data also helps rule out other 
possible origins for the band doubling since the crystal 
structures for Bi2212 and Bi2201 are very similar. For 
example, if the additional band in Bi2212 was a Bi — O 
band or a Cu^d^^^_^^ — 02p hybrid then we would 
expect to see them showing up in Bi2201 as well. The 
lack of the apparent doubling in Bi2201 seems to argue 
against microscopic phase separation, however since the 
mechanism for the phase segregation is not completely 
conclusive, we hold this as a less likely mechanism. 

The /c-space dependence of the splitting of the two 
bands is also consistent with the bilayer splitting effect. 
Symmetry arguments dictate that the coupling between 
the planes should vanish along the (0, 0) (tt, tt) sym- 
metry line and grow as the (7r,0) region is approached, 
exactly as the data in figure 1 shows. More precisely, 
the bilayer splitting has been parameterized by A(fc) = 
0.5*t ±{cos{kxa) — cos{kya))^ which has a maximum split- 
ting of 2tj_ at (tt, 0) [ [l5|Jl6| . Using this equation to pa- 
rameterize the overdoped data of figure 1 and 2, we find 
a, t± oi ^ 55meV and an extrapolated splitting at (tt, 0) 
of llOmeV. This is small compared to the band struc- 
ture prediction of nearly 300meV [ pTf implying impor- 
tant correlation effects present in the system, but it is 
still larger than or comparable to other important en- 
ergy scales {Tc,T*, J) meaning that it is an effect which 
should not be ignored for the physics of these compounds. 

Although we have not yet carried out a full doping 
study, our preliminary data indicates that near-optimal 
or underdoped samples have a smaller splitting effect. 
This smaller splitting, as well as the broader peaks asso- 
ciated with less heavy doping, make the detection of the 
splitting much more difficult, probably explaining why 
this effect has not been observed previously. We also 
note that it would be harder to argue for the existence 
of the splitting from the conventional EDCs because the 
lineshapes are poorly known. For the MDC's which have 
only been used extensively in the past year or so, the 
lineshape is simple and additional peaks can be more 
readily seen, even if the peaks are not as sharp or split- 
ting as large as it is in the data shown here. Addition- 
ally, much effort has been concentrated along cither the 
(0,0) {tt, tt) nodal line where the splitting vanishes, 
or near the (tt, 0) point where the superstructure bands 
contaminate the data and give additional peaks. Hence 
for a clear determination of the bilayer splitting effect it 
is important to study the middle of the zone, as done 
here. 

Also relevant to why it may not have been previ- 
ously detected are the matrix element effects, which may 
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strengthen one of the peaks at the expense of the other. 
For example, for the data of figure 1 and 2, a large out- 
of-plane component of the photon polarization was used, 
which helped to enhance the already-weak B peak. The 
identical data taken with a more in-plane photon polar- 
ization shows a much weaker B peak, as shown in fig- 
ure 3b and 3c, however the fitting result produces same 
amount of the splitting in energy (data not shown here) . 
The major difference between the data with less grazing 
photon incidence angle (case in figure 3) and more graz- 
ing angle (case in figure 1 and 2) is the relative intensity 
of peak A and B which makes the feature B barely visible 
in figure 3b. Of course, the relative strength of A com- 
pared to B can also depend on other parameters, most 
notably the photon energy, for which oscillations in the 
relative magnitude of the branches may exist. 

Although the bilayer splitting scenario is convincing, 
we note that there is increasing evidence for microscopic 
phase separation in BSCCO and related compounds, and 
such an explanation for the band splitting is hard to com- 
pletely rule out However, if this is the explana- 
tion a couple of important restrictions can be made from 
this study: (1) The degree of phase separation should be 
larger for Bi2212 than Bi2201. (2) There would be metal- 
lic regions of two dominant types (possibly with other 
insulating regions) as the MDCs show two main peaks. 
(3) The angular width of each MDC peak is on the or- 
der of 1° HWHM or 0.04^1^^, giving a mean free path 
for scattering I = 1/Ak ~ 25^4 which sets the minimal 
domain size in the sample. 

From our new data, we see that the B piece of FS 
matches the traditional hole-like piece of FS very well 
H-g, while the A_ piece deviates in fc-space significantly 
and approaches M{tt,0). While an extrapolation is re- 
quired, it appears that the A piece probably remains 
centered around Y{t:,tt), i.e. it remains to be hole-like. 
However, it comes close enough to (tt, 0) that a very small 
amount of dispersion of this band could easily push 
it back and forth across (tt, 0), making it cither electron- 
like or hole-like. Experimentally, is varied by changing 
photon energies, which also should affect the intensity ra- 
tio of the two components. Assuming a small amount of 
kz dispersion, this would give a very natural explanation 
for the recent controversy over the FS topology, in which 
the existence of an electron-like portion centered around 
(0, 0) has been suggested as well as denied ||J^. 

Such a kz dispersion would come from a c-axis inter- 
cell coupling, and while this has yet to be directly ob- 
served, wc note that the existence of the intracell cou- 
pling shown here makes the non- vanishing intercell cou- 
pling more plausible. Recent optical measurements on 
Tl-based cuprates do suggest the presence of this inter- 
layer coupling JTot . Similar to the intralcell coupling, 
symmetry arguments tell us that the intercell coupling 



should vanish along the (0,0) — > (7r,0) nodal fine and 
should be largest near (tt, 0), precisely where they would 
need to be to resolve the controversy. However, how to 
connect above argument and Bi2201 FS topology at 33eV 
is still not clear to us. 

In summary, we have shown that in the normal state 
of Bi2212 there is a doubling of the band structure. Near 
the middle of the zone the angular splitting is a little 
more than 2° and the energy splitting is roughly 75meV. 
The /c-dependence of the splitting is consistent with that 
expected from bilayer splitting, and if we parameterize 
this, we obtain a intra-cell coupling t± of ^ 55meV. An- 
other less likely origin for the splitting is the microscopic 
phase separation. Similar results on even more heavily 
overdoped samples have been independently obtained by 
the Stanford group as well [Q. We acknowledge sam- 
ple preparation help from J. Koralek and beamline sup- 
port from H. Hochst, S. Kellar and D.H. Lu. This work 
was supported by the NSF Career-DMR-9985492 and the 
DOE DE-FG03-00ER45809. SSRL and the ALS are op- 
erated by the DOE, Office of Basic Energy Sciences and 
the SRC is supported by the NSF. 
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